Calibrating an interferometric laser frequency stabilization to MHz precision 
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We report on a calibration procedure that enhances the precision of an interferometer based fre- 
quency stabilization by several orders of magnitude. For this purpose the frequency deviations of the 
stabilization are measured precisely by means of a frequency comb. This allows to implement several 
calibration steps that compensate different systematic errors. The resulting frequency deviation is 
shown to be less than 5.7 MHz (rms 1.6 MHz) in the whole wavelength interval 750 — 795 nm. Wide 
tuning of a stabilized laser at this exceptional precision is demonstrated. 
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I. INTRODUCTION 



Fabry-Perot interferometers have advanced to very 
high precision wavelength sensors [HQ- In a Fizeau 
type setup and combined with an absolute frequency ref- 
erence (typically an Helium-Neon laser) they are used 
as precision optical wavelength meters. By means of an 
additional feedback, stabilization and tuning of the laser 
wavelength can be accomplished [3[. However, due to the 
comparatively long signal processing times, the feedback 
bandwidth of these systems is low and limiting the use 
for laser frequency stabilization. 

In this paper we present the optimization and calibra- 
tion of an existing Fizeau type quadrature interferometer 
4] that is used for laser frequency stabilization and allows 
a very wide tuning range of a stabilized laser. Several or- 
ders of magnitude in precision over a wide wavelength 
range of about 45 nm are gained in comparison to an un- 
calibrated setup. The achieved accuracy with an rms de- 
viation of 1.6 MHz is state-of-the-art interferometric laser 
frequency stabilization at high feedback bandwidth. In 
order to characterize and compensate the main observed 
frequency deviations occurring in the wavelength range 
750 — 795 nm we use a fiber laser based frequency comb. 
We take these frequency deviations into account by a cor- 
rection of the wavelength dependence of the refractive 
index up to second polynomial order. This allows for 
future calibration of the setup by employing only three 
interpolation points. For these measurements a Doppler- 
free saturation spectroscopy is easily used as an absolute 



frequency reference. 



II. FUNCTIONAL PRINCIPLE 

Prior to the description of the calibration of the in- 
terferometer, we only briefly review it's functional prin- 
ciple, as it is explained in detail elsewhere Q. Fig.Q] 
shows a sketch of the interferometer setup. A test laser 
beam enters the device through an optical fiber and is 
split into two by means of a wedged beam splitter. The 
beams pass an etalon under a slight relative angle. Two 
pairs of photo diodes, one pair for each beam, are used to 
electronically construct two normalized periodic interfer- 
ometer signals, U x and U y , from the beams transmitted 
and reflected by the etalon [l-Q . The angle between the 
beams results in a tt/2 relative phase shift between the 
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FIG. 1: Patented interferometer setup subject to calibration. 
The interferometer is realized as an etalon of BK7 with a 
geometrical thickness of 50 mm. PD: photo detector, quad: 
quadrature, norm: normalization. See [J] for details. 
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interferometer signals. Due to the etalon's low finesse, 
their frequency dependence can be approximated by sine 
and cosine functions: 



U, 



cos if 
simp 



with ip = 4ttu- 
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(1) 



Hence, U x and U y can be used as quadrature signals de- 
scribing a circular path in the x-y plane. Here, the inter- 
ferometer phase ip depends on the laser frequency v, and 
the thickness L and the refractive index n of the etalon. 
The latter together constitute the etalon's wavelength 
dependent optical path length n (A) • L. Here, en is the 
speed of light in vacuum. We call the frequency interval 
corresponding to a phase change of A(p = 2ir the etalon's 
free spectral range (FSR) , which is approximately 2 GHz 
for the described setup. 

In order to obtain an error signal for frequency 
stabilization of a laser, the interferometer signals 
are compared with electronically generated set signals 
(C4et,x, U S et, v ) = (cos<£ S et,sin</J se t) that depend on a sin- 
gle chosen set phase <p sc t- The difference between tp and 
p sc t is used as a frequency error signal that varies with 
constant slope throughout the free spectral range. This 
allows to use the etalon for frequency stabilization of the 
laser at any point within the FSR. Further, tuning of a 
stabilized laser can be accomplished over a large wave- 
length range including many FSR by alteration of the 
generated p sc t- 

As the interferometer signals are periodic functions of 
the optical frequency, they are ambiguous. Hence, for 
stabilization to a specific target frequency one has to 
predetermine the laser's frequency to an accuracy of bet- 
ter than one half of an FSR. This can be achieved by 
means of a standard wavemeter. Further, tp is subject 
to a drift of n ■ L. We use a reference laser of precisely 
known frequency to keep track of and correct for small 
optical path length changes, that result in an offset of 
the measured phase at the given frequency. In practice, 
as typical interferometer drift rates are on the order of 
100 — 200kHz/min, this measurement has to be carried 
out about every 10 minutes, which can be done automat- 
ically within approximately one to two seconds. With the 
help of these additional measures, each frequency is un- 
ambiguously assigned a phase, and continuous scans over 
an almost arbitrary amount of FSRs can be carried out. 



III. CALIBRATION AND 
CHARACTERIZATION 

To reach an absolute frequency precision on the MHz 
level over a wavelength range of several tens of nanome- 
ters, three calibration steps are performed. 

Firstly, a systematic deviation due to the approxima- 
tion of the exact Airy function describing the interfer- 
ometer signal as a sinusoidal function of the laser fre- 
quency in Eq. (TTJ) has to be considered. This is in practice 




quadrature 
interferometer 



Wavemeter 



FIG. 2: Setup for calibration and characterization of the in- 
terferometer. The frequency comb is a self referenced femto- 
second fiber laser based design (Menlo Systems). A Hydro- 
gen maser and GPS are utilized as time standards. The fiber 
etalon is used to record the data used in the Look-Up-Table 
(LUT). A beat note with a frequency comb is recorded us- 
ing a spectrum analyzer. Application laser: Toptica DLPro, 
Reference laser frequency 384.227981 THz (DFB diode laser, 
TEM Messtechnik), Saturation spectroscopy (Cosy, TEM 
Messtechnik) , Wavemeter: High Finesse WS/7, FI: Faraday 
isolator, S: Shutter, BS: Beam Splitter, PD: PhotoDiode, OF: 
polarization maintaining Optical Fiber, FP: Fiber Port. 



done by calibrating the set phase for a given target wave- 
length throughout one FSR. The calibration is achieved 
by comparing the frequency response of the interferome- 
ter with that of an etalon made of an optical fiber. With 
a fiber length of 2.5 m the FSR of 60 MHz is much smaller 
than the FSR of the interferometer. Therefore, the fiber 
etalon's fringe pattern varies much faster when the laser 
frequency is tuned as compared to the interferometer sig- 
nals U x and U y . This allows to use the fringes of the fiber 
etalon as equidistant markers for frequency intervals. In 
this way the deviations from the linear frequency depen- 
dence of the interferometer phase ip can be measured. 
The maximum nonlinearity found for the quadrature in- 
terferometer is on the order of 2% of one FSR, corre- 
sponding to a maximum of 40 MHz deviation in target 
frequencies within one etalon FSR of 2 GHz. In order 
to compensate this systematic deviation in practice, the 
recorded data are electronically represented in a look-up- 
table (LUT). The LUT is used to modify the set phase 
t^sot for frequency tuning of the stabilized laser. With a 
set phase corrected by use of the LUT we find a remain- 
ing nonlinearity of 0.05% using the same method. This 
corresponds to an accuracy of better than 1 MHz within 
one FSR. It is important to note that this procedure is 
easily repeated and an exact LUT can be created when 
the interferometer is used in a different wavelength range 
at a different refractive index of the etalon. However we 
found that one LUT remained valid when changing the 
frequency of the laser by several hundred FSR. 

In the following we describe a second calibration step 
that enables precise tuning over nanometers throughout 
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a wavelength range as wide as 750 — 795 nm. Over such 
a broad wavelength range the dispersion of the etalon's 
medium is significant. The index of refraction can be 
described by the Sellmeier equation Q with the specific 
material parameters for the BK7 ctalon used in this case. 
However, in order to be able to calculate tp se t based on 
Eq.rfjwith a precision corresponding to 1MHz, n(X) ■ L 
has to be determined with a relative error on the or- 
der of 10 -9 . In this calibration step we obtain a pre- 
cise measurement of the optical path n(X) ■ L with the 
measurement setup shown in Fig.[2J For an absolute fre- 
quency reference we use a reference laser at the D2 op- 
tical transition line of 87 Rb at 384.227981 THz gj]. This 
laser is stabilized to the cross-over signal between the 
5 2 Si/ 2 , F = 2 — > 5 2 P 3 / 2 , F = 2,3 transitions by means 
of a saturation spectroscopy. For a precise determina- 
tion of the frequency of the application laser controlled 
by the interferometer setup we employ a frequency comb. 
The frequency is inferred by recording the beat note of 
the application laser with the comb, in combination with 
wavelength measurement carried out by a high resolution 
wavemeter. In order to collect data over a wide wave- 
length range the following steps have to be repeatedly 
performed: 

• The application laser is frequency stabilized to the 
interferometer with a given phase ip sc t at a given 
wavelength. 

• The fiber etalon is used to create an LUT as ex- 
plained in the first step. 

• The reference laser is used to perform a measure- 
ment of the interferometer phase offset. 

• A precise measurement of the application laser fre- 
quency is carried out by means of the frequency 
comb for various ip se t around the given wavelength. 

The results of this procedure for different wavelengths are 
summarized in Figs. OJa) and|4l To obtain the remaining 
frequency deviation after applying the previous calibra- 
tion step we process the data in the following way: We 
first use Eq.[T]to convert the measured frequencies to in- 
terferometer phases and compare these to the set phases 
(p se t. This resulting wavelength dependent phase devia- 
tion can be expressed as a frequency deviation /d ev by 
scaling with the etalon's FSR. The resulting data exhibit 
systematic deviations in the frequency control of the ap- 
plication laser on the order of several ten MHz over a 
wavelength range of 45 nm. This corresponds to a rela- 
tive error on the order of 10~ 7 . 

In order to further characterize the performance of the 
device we show in Fig. [4] the same data resolved for the 
two wavelength ranges that are indicated in Fig.|3](a): 

Fig.|4](a) shows the frequency deviation occurring 
within one FSR. This gives a measure on how well we 
are able to calibrate the relation between interferometer 
phase and frequency by means of the LUT. For a sin- 
gle scan we find small rms deviations. However, for re- 
peated calibrations and subsequent scans over one FSR 



the precision of the device is limited by the reproducibil- 
ity to an rms value of 0.9 MHz. This can be attributed 
to frequency offsets consistent with the linewidth of the 
reference laser. 

Fig. 2](b) demonstrates the capability of precise fre- 
quency tuning over a range as broad as of 225 FSRs. For 
this measurement a distributed feedback diode is used as 
an application laser, which was continuously stabilized 
throughout the scan. Again using the dataset shown 
as squares for calibration, we find an rms deviation of 
0.7 MHz for the dataset represented with circles. Thus, 
the calibration step described above proofs valid for the 
complete wavelength range over approximately 1 nm. 

To calculate the interferometer phase from the mea- 
sured frequency in the above mentioned analysis we use 
the refractive index as described by the Sellmeier equa- 
tion with coefficients specified for the BK7 in use, to- 
gether with an estimated etalon length L. The optical 
path length has been estimated such that the observed 
deviations are minimal throughout the full wavelength 
range shown in Fig. 0(a). As a constraint the optical 
path length was chosen such that the phase offset mea- 
sured at the wavelength of the reference laser was repro- 
duced. With this the second calibration step, in which 
n{X)-L is determined to a very high precision, is achieved. 
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FIG. 3: (a) measured frequency deviation fdev over the whole 
accessible wavelength range. Circles and squares are separate 
datasets taken on two consecutive days. One second order 
polynomial fit to both data sets is shown. The measurements 
for the frequency ranges 1 and 2 are shown enlarged in Fig. [4] 
(b) residual of a fit to both datasets shown in (a). A vertical 
line marks the reference laser frequency in both graphs. 
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FIG. 4: Frequency ranges 1 and 2 shown with higher resolu- 
tion, (a): Frequency range 2. Three sets of data points at 
frequency 396 THz within one FSR of the etalon are shown. 
Circles and squares are taken using different offset phase mea- 
surements but the same LUT. For the triangles a new LUT 
was generated and the offset phase was newly measured. A 
linear drift of f30kHz/min obtained from repetitive measure- 
ments at the reference wavelength has been subtracted from 
all datasets for clearer visibility of the frequency-dependent 
error, (b): Frequency range 1. Two datasets close to the 
reference laser frequency (vertical line) are shown. To record 
this data, the reference laser (DFB diode laser) was controlled 
by the quadrature interferometer and interchanged with the 
application laser (Toptica DLPro). 



IV. DISCUSSION 

It can be seen in the data presented in Fig.[3](a), that 
the frequency deviation is found to be dominantly of sec- 
ond order. This could occur due to a second order er- 
ror in the wavelength dependence of the refractive in- 
dex described with the specified Sellmeier parameters. 
In the second calibration step, we have chosen an op- 
tical path length which minimizes the overall frequency 
deviation in the whole wavelength range and leads to 
zero frequency deviation at the reference laser frequency. 
With this choice of n ■ L, we compensate the zeroth or- 
der error in the description of the refractive index and 
obtain corrections in higher orders to the wavelength de- 
pendence of the refractive index. We find relative er- 
rors of 4.4 x 10~ 5 ± 1.6 x 10~ 6 in the first order and 
8.0 x 10~ 2 ± 2.9 x 10~ 3 in the second order. These de- 
viations are consistent with a relative uncertainty of the 
refractive index described by the Sellmeier parameters of 
2 x 1CT 7 . 

As we have shown that a correction of the refractive 
index up to second order serves very well to calibrate the 
setup, it will from now on only be necessary to measure 
the deviation at three roughly equally spaced frequencies 
within the wavelength range. The calibration is then ob- 
tained by using a second order polynomial fit to these 
data points in the way as described above. In the exam- 
ined wavelength interval, Doppler-free spectroscopies of 
optical transitions in rubidium and potassium are readily 
available and can serve as absolute frequency references. 
This eliminates the need for a frequency comb to carry 
out the last calibration step. In addition, employing data 
at just three wavelengths allows for efficient recalibration 
in order to compensate possible slow systematic drifts on 
the time scale of days. 



The characterization obtained with the data shown in 
Fig-El(a) can be utilized for a third and last calibration 
step. For practical purposes we use a best fit to the data 
by a second order polynomial that serves as a calibration 
curve and is shown in Fig. 0(a). The residuals of this 
fit are shown in Fig.[3](b). Based on the obtained curve 
we can implement the calibration of the interferometer 
setup by applying a wavelength dependent correction for 
the refractive index up to second order. This effectively 
changes the etalon's optical path length used in our calcu- 
lation of the interferometer phase for a target frequency. 

To further demonstrate the effectiveness of this cali- 
bration we use a different second order polynomial fit 
function that has been obtained from the dataset repre- 
sented as circles in Fig.[3](a). In this manner we find a 
maximum frequency deviation of 5.7 MHz (rms 1.6 MHz) 
in the dataset represented as squares for the fully cali- 
brated interferometer. Comparison of two datasets ob- 
tained on consecutive days (circles, squares in Fig. [3]) in- 
dicate a small systematic drift that is not further inves- 
tigated within this work. 



V. CONCLUSION 

In this paper we describe a calibration procedure for an 
interferometric laser frequency stabilization in the wave- 
length range 750 — 795 nm. Several calibration steps 
are taken resulting in an exceptional precision with an 
rms deviation of 1.6 MHz. The maximum observed de- 
viation of 5.7 MHz corresponds to a relative accuracy 
of 1.4 x 10~ 8 . While the accuracy within one FSR is 
increased by approximately one order of magnitude by 
the first calibration step, a calibration over a wide fre- 
quency range allows the accurate determination of the 
interferometer's optical path length. Only this increases 
the precision of the device from hundreds of MHz to the 
MHz level over the wide frequency range. The calibra- 
tion method presented here is expected to be applicable 
throughout the 250 nm wavelength range of the interfer- 
ometer setup, which is limited by the choice of the single 
mode optical fibers used in the system. Further, contin- 
uous high bandwidth laser frequency stabilization over 
a wide tuning range is supported. This makes the pre- 
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sented laser frequency stabilization an optimal choice for 
the spectroscopy of molecules over a large wavelength 
range. 

The accuracy of the device is enhanced, if frequency 
stabilization in a comparatively smaller wavelength range 
of approximately 1 nm is required, as shown in Fig.|4](b). 
If in this reference range a reference laser is available, the 
third calibration step can be omitted. If for the respec- 
tive wavelength range no laser spectroscopy is available 
as a reference, a beat note with a known frequency comb 
mode can be used as a reference. This allows then to con- 
tinuously frequency stabilize the application laser while 
tuning over a range spanning hundreds of etalon FSRs 



of 2 GHz. In comparison, a direct beat note of the ap- 
plication laser with the comb mode can be applied con- 
tinuously only throughout a fraction of the narrow mode 
spacing of the comb. 
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